Planar Sheets Meet Negative Curvature Liquid Interfaces 
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If an inextensible thin sheet is adhered to a substrate with a negative Gaussian curvature it 
will experience stress due to geometric frustration. We analyze the consequences of such geometric 
frustration using analytic arguments and numerical simulations. Both concentric wrinkles and eye- 
like folds are shown to be compatible with negative curvatures. Which pattern will be realized 
depends on the curvature of the substrate. We discuss both types of folding patterns and determine 
the phase diagram governing their appearance. 



Geometric frustration occurs in wrapping a spherical 
Mozartkugel ("Mozart sphere") with planar foil [H-Q- 
The extra circumference at the edge of the planar sheet 
compared to the spherical substrate to which it is con- 
forming gives rise to ridges, narrow deformed regions that 
occupy a small fraction of the total available volume and 
along which the energy is focused Q]. Recently the frus- 
tration of a thin circular elastic sheet of ~ 1mm size 
covering the cap of a spherical droplet has been studied 
At first fine radial wrinkles appear at the edge of 
the sheet and then become unstable to localized folds as 
the size of the droplet decreases. Once expects the sit- 
uation to be quite different for a sheet conforming to a 
negative curvatures surface - in this case the edge of the 
sheet is stretched tangentially as opposed to being com- 
pressed on spherical geometry. A completely different 
frustration pattern on the planar sheet is thus expected. 
In this paper, we study the wrinkle/fold structure in an 
elasto-capillary system on a flat sheet conforming to a 
negatively curved substrate. 

The system we treat has two parts: an inextensible 
thin elastic hydrophilic sheet and a saddle-like fluid in- 
terface with negative curvature. The size of the clastic 
sheet is t aken to be much bigger than the elasto-capillary 
length y/ n/cr, so that surface tension a dominates over 
the bending rigidity k 0. Note that a standard sheet 
of paper can be regarded as inextensible and be used to 
demonstrate the bending of a thin elastic sheet Q • When 
such an inextensible elastic sheet is placed on a negative 
curvature liquid interface the capillary force pulls the pla- 
nar sheet into full contact with the liquid interface. 

Complete contact between a planar sheet and a curved 
fluid substrate by capillary forces introduces a wrin- 
kle/fold pattern, which redefines the metric of the pla- 
nar sheet according to the curvature of the background 
geometry. The modification of the metric leads to the 
change of shape. This concept has been utilized to de- 
sign responsive buckled surfaces [l^ . The curved sur- 
face endows its metric to the flat sheet via their full con- 
tact. The inherited metric on the flat sheet as well as its 
elasticity determines the wrinkle/fold structure. Which 
wrinkle/fold patterns are compatible with negative cur- 
vatures? The art of origami provides some inspiration. 
It is shown in origami that regular concentric or square 
wrinkles with alternating peaks and valleys on a piece of 
paper can induce a negative curvature, as demonstrated 




FIG. 1. (color online) Wrinkles of alternating peaks and val- 
leys on a piece of paper redefine the metric, automatically 
bringing the flat paper to a saddle-like shape. It is shown in 
the text that the effect of concentric wrinkles is equivalent to 
inserting a wedge of some angle in the sense of redefining the 
metric, (a) is excerpted from Ref. @| , and (b) is from Ref. [sj 
where the method to make it is described in detail. 



in Fig. [Q 

In what follows we will prove that the effect of concen- 
tric wrinkles is equivalent to inserting an angular wedge 
in the sense of redefining the metric. The wavelength 
and amplitude of the concentric wrinkles determine the 
angle of the wedge. The wrinkled sheet is parameter- 
ized as x(r,9) = {rcos#, rsin0,afccos(fcr)}, where 
is the amplitude of the wrinkles and k is the wavenum- 
ber k = 2ir/\. The nonzero components of the metric 
tensor are g\\ — 1 + x 2 sin 2 (fcr) and (722 = r 2 , where 
x = a^k. The imposed wrinkles transform the original 
wrinkled shell into a new surface denoted by £ that we 
take coincident with the curved substrate. On this new 
surface the metric is redefined such that the geodesic dis- 
tance (denoted as r) from the center of the disk to the 
first peak is tq = A. The corresponding real distance 
on the original sheet is Iq = j Q X ^/gudr = \E (x), where 
E (x) = f n 2w dy\/ 1 + x 2 sin 2 y. As — >• 0, In — > A, as ex- 



/„" dyy/l + x 2 sir 
pected. Due to the inextensibility of the paper model, the 
mapping from the originally flat sheet to a wrinkled shape 
is isometric and thus length-preserving. The perimeter 
of the circle with radius r in the new surface S is thus 
C (r) — 2irl (r) = 2nr+ra, where Z(r) = l^kr j (2tv). Here 
a is interpreted as the angle of the inserted wedge: 
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The positive sign of a indicates that concentric wrinkles 
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are equivalent to inserting a wedge in redefining the met- 
ric. The inserted wedge buckles a flat disk to a saddle- 
like shape with negative curvature in three-dimensional 
Euclidean space [111 ], as do concentric wrinkles. These 
two ways of introducing negative curvature - either in- 
serting material or imposing concentric wrinkles - are 
related via the expression for C(r). The first method 
changes the perimeter without changing the radius, while 
it is the opposite for the second method. By expand- 
ing the expression for l in terms of small x, one finds 
C{r) = 2irr + ^k 2 x 2 r 3 + 0{x). By inserting this into 
K G = lim r ^o 3 [2?rr - C (r)] / (nr 3 ) PJ], the curvature 
at the center of the surface £ is found to be 



Kg = 



,2.4 



fc . 
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Eq. ([2]) shows that increasing the amplitude or the fre- 
quency of the concentric wrinkles results in greater cur- 
vature of the surface, with greater sensitivity to the fre- 
quency. The curvature of the background geometry de- 
termines the amplitude and wavelength of the concentric 
wrinkles according to the product a 2 k 4 . 

It is worth mentioning that by increasing ( "growing" ) 
the radius of a flat disk while keeping the perimeter 
invariant, we get a positive curvature surface, which 
may buckle to various patterns depending on its elas- 
ticity. This elasticity paradigm has been used to explain 
the phyllotactic patterns of Fibonacci-like sequences on 
plants [l3j . Wrinkles of concentric squares can buckle a 
square piece of paper to a beautiful hyperbolic parabola 
with negative curvature as shown in Fig.fljb). Prescribed 
metrics via the design of wrinkles can even transform a 
flat piece of paper into a rich variety of structures, in- 
cluding the DNA double helix fill ]. 




FIG. 2. (color online) A branch-like fold pattern appears on a 
flat paper disk by wrapping it on a negative curvature surface 
(left figures). The folds in the red rectangles are equivalent 
to removing eye-like areas as shown in the upper right figure. 
The buckled shape (the lower right figure; the leaf-like ob- 
ject therein is the removed material) of a flat disk due to an 
isolated fold is obtained by "closing the eye". 

In addition to wrinkles, localized folds can also change 
the metric [HI]. By attaching a paper disk to a negative 
curvature surface one finds typical branch-like fold pat- 
terns as shown in Fig. [5] (the left two photos are the same 



deformed paper disk from different perspectives). A light 
beam directly illuminates the paper disk from above, so 
the folds are seen clearly as black curves. The folds on 
the sheet can be roughly classified as principal ones (in 
red rectangles) and fine ones (barely seen above the red 
rectangle in the left lower figure). Their role in "screen- 
ing" curvature is similar to top olog ical defects in crys- 



talline order on a curved surface [161 ] . The feature of the 
observed folds is that the amount of the folded material 
decreases towards their ends. Similar folds are also found 
in the interior side of a bent tube where the curvature is 
negative [l?}- We analyze an isolated fold to illustrate 
that it is compatible with negative curvature geometry, 
which is analogous to the compatibility of a seven-fold 
disclination with negative curvature geometry. The ef- 
fect of such folds is to remove an eye-like area from a flat 
disk. The buckled shape (see the lower right figure in 
Fig. [2]) due to the fold, or equivalently the removal of an 
eye-like area, is obtained by "closing the eye" . The cur- 
vature of the buckled shape is negative; the disk curves 
up along the fold and curves down along the orthogonal 
direction. The removal of an eye-like area is equivalent 
to inserting a wedge, because more material is removed 
at the center than on the edge. The profile of the eye- 
like fold can be determined by the curvature of the back- 
ground geometry. The perimeter of a circle with geodesic 
radius r on a disk with an eye-like removed is esti- 
mated as C(r) ss 27rr + 4 (h (0) - h (r)), where 2h (r) is 
the width of a fold. For a small-sized fold (in compar- 
ison with yJl/\K G \), one finds h(r) = h(0) + j^Kcr 3 
by inserting C (r) into the expression for the Gaussian 
curvature. From h(r = L/2) = 0, we see that the length 
of the fold is controlled by the curvature according to 

1 /3 

L = 2 ^ ^Kg] ) ' ^ n con trast, it is interesting to note 
that reversed eye-like folds are found in one's palms. It 
seems that the two main lines - head line and heart line 
following the terminology of palmistry - are compatible 
with positive curvatures as the width of these lines in- 
creases from the center to the edge of a palm. 

We now examine the transition between wrinkle and 
fold patterns. In the regime of large surface tension, 
the ground state of an elasto-capillary system is dom- 
inated by the surface energy difference before and af- 
ter a planar sheet is attached to the liquid interface: 

Ai* 1 = CTLA^sheet + &LS -^coverage ~ O'LA^substrate = 

— (TLsAcoverage + const. (T/j is the surface tension be- 
tween phases / and J with L, A, and S standing for 
liquid, air, and elastic sheet, respectively. A coverage is 
the area of the liquid substrate occupied by the clastic 
sheet, which is smaller than the area of the sheet A s h ee t 
due to its deformation. A su b stra te is the sum of the occu- 
pied and unoccupied substrate areas, which is a constant. 
The surface energy turns out to depend only on A coverage ; 
the larger it is, the smaller the energy is. Unlike a planar 
sheet on a positive curvature surface, the deformation at 
the edge of a planar sheet on a negative curvature surface 
may be ignored. Therefore, the optimal contour shape of 
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FIG. 3. (color online) The phase diagram of a deformed sheet 
on a negative curvature surface, z (x, y) = ax' 2 — f3y 2 . Wrin- 
kles (squares) on the sheet occur near the isotropic region 
(a//3 ~ 1) while folds (triangles) are found in the highly 
anisotropic regions (a//3 far from unity). For very small a 
and ft the sheet is compliant with the adhering surface. This 
region is designated as smooth (circles). Filled symbols indi- 
cate parameters for which simulations were performed. The 
phase diagram is symmetric with respect to the a — f3 line. 
Mirror images of the simulated points are represented as open 
symbols. 



Here a,/3£ {x,y}, E is the three-dimensional Young's 
modulus and Y = Et with t being the sheet thick- 
ness; u a fj = ^ (g a fs — Hap) is the strain tensor with 
Vapidap) being the metric tensor of the reference (de- 
formed) state. The reference metric is assumed to be 
flat, i.e., g a p = 8 a p. b a p is the second fundamental form 

[20l ]. dA = y^ldxdy is the area element and v is the 



Poisson's ratio. Finally, 



J , Q7 and bP = g Pl b ai . 
With the mean curvature H = 1 6^ and the Gaussian cur- 
vature K G = det E b = f dAn (2H 2 - (1 - v) K G ), 
where the bending rigidity is k = {i-u' 1 ) ■ Since the 
deformations of the sheet do not change its topology, the 
Gauss-Bonnet theorem [2(| ensures that J dAK G is con- 
stant and the Gaussian curvature term can be omitted. 

We find the relaxed shapes of the sheet by perform- 
ing simulated annealing Monte Carlo (MC) simulations 
of a discrete triangular mesh. The discrete form of the 
stretching energy Eq. © is [HI, S3] 
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where F = g g — /, At is a mesh triangle area, and 
the sum is carried out over all triangles, g and g are the 
discrete counterparts of the reference and actual metric 
tensors, respectively, whose elements are the scalar prod- 
ucts of the two vectors spanning each triangle before and 
after the deformation. By assuming that the sheet has 
no spontaneous curvature, the bending energy (Eq. Q) 
of the discrete mesh is [23, r24j 



a deformed sheet on the liquid substrate is the one that 
maximizes the coverage area while keeping the perimeter 
fixed. This is exactly the classical isoperimetric problem, 
in this case on a negative curvature surface. For constant 
curvature surfaces, the classical isoperimetric solution in 
the Euclidean plane is also valid with the circle in E 2 
being replaced by a geodesic circle [lU. On a general 
surface with varying negative curvature, there is no ex- 
act mathematical result available. The physical picture, 
however, is rather interesting: a deformed sheet fully at- 
tached to a curved liquid substrate will migrate to the 
region where it can extend as far as possible to maximize 
the contact area; the driving forces are the capillary force 
and the release of the bending energy in this curvature- 
driven migration process. 

In order to further explore the deformation patterns 
we have performed numerical simulations of a thin clastic 
sheet adhered to a substrate with negative curvature. For 
an isotropic material the stretching and bending energies 
are given as [l9| 



E b = — dA 
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where V r; is the gradient with respect to the position 
of vertex i and the sum is over all vertices. Aj = 
Ax is the vertex area, with £7^ being the ver- 
tex "star" , i.e. the set of all triangles that share vertex 
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and are the coordinate vectors of the vertices of T. 
We omitted (T) in the superscript of vector since it 
is shared by all T E fij. N; = V r<i Vi is a volume gradi- 



(rf ) - nj x (rf) - r<) 



where r,-, r 
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ent, with the volume |25| V* 



£t(« 



(T) 



The sheet is assumed to adhere to the substrate via a 
potential, E a dh. = \l~Yliid 2 Ai, where di is the shortest 
Euclidean distance between the vertex i and the adhering 
surface. 

The surface mesh with pa 2 x 10 4 triangles wasgen- 
erated by constructing a Delaunay triangulation [26[ of 
« 10 randomly but evenly distributed points on a disk 
of radius R ~ lOOao, where eto is the average distance 
between two neighboring points. The initial state for 
each simulation was constructed by deforming the pla- 
nar mesh to perfectly comply with a prescribed shape 
of the substrate, modeled as a hyperbolic paraboloid and 
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FIG. 4. (color online) Snapshots of the relaxed conformation 
of a sheet adhered to a negative curvature substrate mod- 
eled as a hyperbolic paraboloid surface, z (x,y) = ax 2 — f3y 2 . 
For a substrate with a highly anisotropic shape (a = O.la^" 1 , 
/3 = O.Sa^ 1 ) one observes folds (a); sheets adhered to an 
isotripic substrate (a — /3 = 0.5a,g 1 ) develop wrinkles (b). 
Mean curvature, H, is measured in units of Oq . 



parametrized as z (x, y) = ax 2 — f3y 2 . A deformation pat- 
tern typically emerged within the first 5 • 10 4 MC sweeps 
and was further relaxed for an additional 2.5 x 10 s sweeps. 
In all simulations the energy scale is set by the bending 
rigidity k and we set Y — 10 3 ft/a§ (corresponding to the 
thickness t w 0.03ao), 7 = IOac/oq, and v = 1/3. 

In Fig. [3] we show a phase diagram of the wrin- 
kle and fold patterns for a, (3 <E [0, ...,0.5], measured 
in units of a^ 1 . We note that the Gaussian curva- 
ture of a hyperbolic paraboloid surface is Kq (x, y) — 
—Aaptj (l + Ax 2 a 2 + 4y 2 {3 2 ^J and has a maximally neg- 
ative K G lax = —4af3 at x = y = 0. For small values 
of a, P < 0.2a^ the Gaussian curvature is very small, 
I .Kg I £ 0.15a^" and the strong adhesion forces prevent 
the sheet from deforming; instead it remains smooth and 
perfectly compliant to the substrate. As a and /3 in- 
crease two distinct deformed states form, folds and wrin- 
kles. These can be distinguished by the isometry con- 
dition as follows. The inextensibility of a sheet sets an 
upper limit 2_R, the diameter of the planar sheet, as the 
maximal geodesic distance between two arbitrary points 
on a sheet contour. If the determined sheet contour has 
the maximum geodesic diameter equal to (smaller than) 



2R, then the sheet is recognized as having folds (wrin- 
kles). An isotropic Gaussian curvature (a = /?, since 
Kq = f {a 2 x 2 + (3 2 y 2 ) [Hf) imposes either isotropic tan- 
gential stretching or isotropic radial compression on the 
sheet, which is expected to result in wrinkles. As the 
Gaussian curvature grows increasingly anisotropic (the 
ratio a /P deviates from unity), the imposed anisotropic 
stretching and compression on the sheet are expected to 
generate folds which are themselves anisotropic objects. 
This shows that wrinkles occur near the isotropic region 
while folds arise in the highly anisotropic regions. A typ- 
ical configuration with folds is shown in Fig. HJi, while a 
typical wrinkle pattern is shown in Fig. |4)d. Finally, in 
the highly anisotropic case, a <C /3, Kq ~ and the sub- 
strate is nearly cylindrical; in this case the deformation 
of the sheet is nearly isometric and no wrinkles or folds 
form. 

In conclusion, using geometric arguments and numer- 
ical simulations of a thin-sheet elastic model we study 
the curvature-driven wrinklc/fold patterns of a flat sheet 
adhered to a negative curvature substrate. We analyze 
two types of structure, concentric wrinkles and eye-like 
folds, that are compatible with negative curvature liquid 
substrates and discuss the transition between these two 
states driven by the anisotropy of the background geom- 
etry. The ability of wrinkles / folds to deform a plane to a 
curved surface may find potential applications. Consider 
a flat sheet with a pre-designed wrinkle/fold pattern like 
the lines on a palm. By controlling the on/off status of 
the wrinkles/folds, a planar sheet can be programmed to 
buckle to a desired shape. This may lead to potential 
applications in maximizing sunshine harvest bydesign- 
ing the shape of ultra thin flexible solar cells [27]]. Our 
study also sheds light on the reverse problem of attaching 
a curved shell to a flat substrate, e.g., the adhesion of a 
cell on a flat substrate Q . 
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